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We present new results for the Kondo lattice model of strongly correlated electrons, in 1-, 2-, and 
3-dimensions, obtained from high-order linked-cluster series expansions. Results are given for varies 
ground state properties at half-filling, and for spin and charge excitations. The existence and nature 
of the predicted quantum phase transition are explored. 

I. INTRODUCTION 

The Kondo lattice model, described by the usual Hamiltonian 

H = -tY,{clc,, + h.c.) + jY,S,-s, (1) 

{ij)a i 

represents a band of conduction electrons, interacting via a spin-exchange term with a set of immobile s = 5 spins Si 
(f electrons). l-. 

The model has been extensively studied in connection- with a class of materials known as "Kondo insulators"Ej 
(J > 0), and in connection with the manganites (J < 0)0. Despite the apparent simplicity of the model, in which 
neither the conduction electrons nor localized spins interact directly among themselves, the spin exchange leads to 
a strongly-correlated many-body system. No exact results are known for cither ground state or thermodynamic 
properties for general J /t, in any dimension. 

The model incorporates two competing physical processes. In the strong-coupling (large | J|) limit, the conduction 
electrons are "frozen out" via the formation of local singlets (J > 0) or triplets {J < 0). In either case there will be a 
gap to spin excitations and spin correlations will be short ranged. On the other hand, at weak coupling, the conduction 
electrons can induce the usual RKKY interaction between localized spins, giving rise to possible magnetically ordered 
phases with no spin gap and long-range correlations. In one dimension there will be a smooth crossover from large 
I J| to small I J| behaviour, but in higher dimension a quantum phase transition is expected. 

A great deal of work has been carried out on thej-one dimensional model, using a variety of analytic and numcydjca.1 
methods, and we refer the reader ta a recent reviewO. In higher dimension there have been mean-field approachesOuEl, 
quantum Monte Carlo calculationsQ, and a series expansion studyB. These studies, which are all for the half-filled case, 
conclude that a quantum phase transition, at which the spin gap vanishes continuously, occurs at {J/t)c — 1.45 ±0.05 
in the 2D J > case, while Refs. ||J| give {J/t)c — 1.833,2.0 respectively for the 3D J > case. There have not 
been, to our knowledge, any similar studies for the case of ferromagnetic coupling. 

Our aim in this paper is to study the Kondo lattice model in 1, 2 and 3-dimensions via series expansion methods. 
We have considerably extended the calculations of Ref. ||, by obtaining longer series, by using also expansions about 
the Ising limit, and by studying also the energies of elementary excitations. 

Linked-cluster series ejxpansions have been used successfully for many years to study strongly interacting lattice 
models. A recent reviewO describes the basic approach and some of the results which have been obtained. The method 
is applicable in any dimension, is particularly suited to locating critical points and is free from finite size corrections 
or minus sign problems which hamper other numerical approaches. On the other hand good convergence may be 
limited to particular regions of the phase diagram. 

The Hamiltonian is written in the generic form H = Hq + XV where Hq has a simple known ground state. The 
remaining term(s) in H are treated perturbatively, to high order. In this way the ground state eneiigy, correlations, 
susceptibilities, etc, are expressedpas-power series in A. These are then analysed by standard methodsE]. An extension 
of the basic linked-cluster methodBllil allows the computation of the full dispersion relation for elementary excitations, 
which can yield energy gaps. 

For the present model the simplest choice is to take Hq — JJ2i Si ■ s^, a sum of single-site exchange terms. The 
unperturbed ground state is then a simple product state of dimer states. This is the approached used in Ref. ^, and 
also our first method here. We refer to these as "dimer expansions" . An alternative is to write the exchange term as 

jT}^:^! + liS+sr+Sr4)] (2) 

i 

and to take only the first term as Hq, In such an "Ising expansion" both the spin- fluctuation and hopping terms are 
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treated perturbatively. To remove the degeneracy in Hq, we add following two terms for conduction spins into Hq 

J'Y.{sls^ + 1/4) + hY,[{-irst + 1/2] (3) 

and subtract them from the perturbation term, so the overall Hamiltonian is 

H = Ho + XV 

Ho = j5](5fsn + J'E(sf4 + l/4) + /i5][(-l)X + l/2] 

i (ij) i 

i (ij) i {ij)a 

Series in power of A are computed for given values of J, i, J' and h, and extrapolated to A = 1 where the original 
Hamiltonian is recovered. Such expansions are appropriate for magnetically ordered phases, although they can also 
yield accurate results in other cases. 

In the remainder of this paper we will present and discuss our results for the ID case (Section II), for the 2D square 
lattice and for the 3D simple cubic lattice (Section III). An overall summary is given at the end. 



(4) 
(5) 



II. THE ID KONDO LATTICE MODEL 



Using the dimer expansion approach we have computed series for the ground state energy in the form 

oo 

Eo/NJ ^Y.^s{tlJY (7) 

s=0 

and the coefficients, to order 20, are given in Appendix A. Our coefficients agree exactly_with Ref. ||, and add three 
new terms (odd coefficients vanish for this series). Integrated differential approximantal3 are used to evaluated the 
series for particular t/J, and the resulting energy is shown in Figure ^ The different approximants agree well up 
to t/ J ~ 1J2, but then splay outwards. We also show, for comparison, the energy obtained from an early DMRG 
calculationli3. 

Next we turn to the spin excitations. In the strong-coupling limit a spin excitation corresponds to a spin-triplet at 
one site, which is able to propagate coherently via the conduction electron hoping term. The dispersion relation can 
be expressed in the form 

As(fc) = ^i„(A)cosnfc (8) 

n 
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FIG. 2: The triplet spin excitation spectrum (a) and one-hole (quasiparticle) excitation spectrum (b) for ID case. 
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where the quantities tn(A) are expressed as power series in A = t/J. For the ID case we have computed these up to 
n = 9 (order 18 in A), and for the interested reader we provide this data in Appendix A. Figure |^(a) shows the triplet 
spin-excitation energy vs k, for value oit/J = 0.25,0.4,0.5. For t/J = the excitation will, of course, have energy 
J and will be dispersionless. Increasing the hopping amplitude gives increasing bandwidth, with the energy at A: = 
raised slightly and a minimum at fc = tt. We are unaware of any previous reported calculations of this dispersion 
relation, apart from the second order result given in Ref. |^. 

From eq. ^ at fc = tt we obtain a series in A for the spin gap, which is again evaluated using integrated differential 
approximants. Results are shown in^^igure I The series is well rConverged up to t/J ~ 1.1. For comparison we 
show spin-gaps calculated by DMRGOj and a mean-field approach^. Agreement with DMRG is excellent over the 
range shown, while the mean-field method appear to seriously under estimate the size of the gap. All of the results, 
including ours, are consistent with a spin-gap which decreases rapidly but does not vanish until J = 0. 

Next we consider the so-called "quasiparticle" excitation, which we prefer to call a one-hole excitation. This 
corresponds to the removal of an electron from the half-filled band and thus, in the strong-coupling limit, to a single 
localized spin on one site with singlets on the others. For t = the energy gap is thus 3 J/4. For the ID case we have 
computed this series up to order 12 in A. Figure ||(b) shows the one-hole excitation spectrum for values t/ J = 0.25, 




2.5 



FIG. 3: The triplet spin gap As/ J and quasipartiete gap Aqp/J versus t/J obtained from different integrated differential 
approximants. The points are the results of DMRGEj for As/ J (points with errorbars) and Ac/2 J (solid circles). The short 
dashed line is the result of a mean-field approachtfl. The inset gives a logarithmic plot for As/t. 
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0.4 and 0.5. The minimum occurs at fc = and the bandwidth seems roughly proportional to t/ J . We are, again, not 
aware of any previous calculations of this dispersion curve. The series at fc = allow us to compute the one-hole gap, 
and our results are plotted in Figure |[ A notable feature is that A^p becomes approximately constant for t/ J > 0.5. 
We know of no previous calculations for the one-hole gap apart, again, from the second order strong-coupling result in 
Ref. ^. There is yet another gap, the "charge gap" , which corresponds to an excitation in which the system remains 
half-filled, but with a doubly occupied site and an empty site. We, are not able to compute this via series, at this 
stage. However, in the strong coupling limit Ac = 2Agp (this is valia3 to at least second-order in t/J). The charge gap 
has been computed by DMRGE3, and we show in Fig. |3[ the result for Ac/2. Evidently for larger hopping parameter 

AJAqp < 2. 

III. THE SQUARE AND SIMPLE CUBIC LATTICES 

In two or more spatial dimensions it is believed that the Kondo lattice model has a true quantum phase transition 
at some (J/t)c, between a gapped spin-liquid uhase and a small J magnetically ordered gapless phase. For the 
square lattice quantum Monte Carlo simulationa3 provide strong indications of a transition at (i/J)c — 0.69, while a 
bond-operator mean-jfield theorjfl gives {t/J)c = 0.664, 0.546 for the square and simple cubic lattices, respectively. A 
previous series studyO has given (t/ J)c — 0.7, 0.5 respectively. These latter estimates are relatively imprecise, and it 
seemed worthwhile to investigate this using longer series. 

We have derived dimer series for the square lattice for the following quantities: ground state energy Eq, antiferro- 
magnetic spin susceptibilities for both local and conduction spins {xi and Xc), and the triplet spin excitation spectrum 
(all to order 12), and the one-hole ( "quasiparticle" ) excitation spectrum (to order 11). This adds two non-zero terms 
to the results of Ref. ||. The excitation series are new. In addition we have computed, for the first time, Ising 
expansions for the ground state energy and for the staggered magnetizations (for both local and itinerant spins) to 
order 13. Series, to the same order, have been derived for the simple cubic lattice for all of the same quantities, except 
for the excitations, where we have only computed the minimum gap rather than the full spectrum. The dimer series 
for ground state energy Eq, antiferromagnetic spin susceptibilities for both local and conduction spins {xi and Xc), 
and the minimum triplet spin gap are given in Table I. Our results agree completely with those of Ref. ^ for the 
square lattice but disagree for the simple cubic lattice susceptibilities beyond the 4th term. We are unable to resolve 
this, but we confident that our results are correct. 

We first show, in Figure ^, our estimate of the ground state energy, as a function of t, obtained from both dimer 
and Ising expansions. Both series converge well for small t, but the Ising expansion has better convergence for larger 
t. There are no previously reported data for Eq. There is no evidence of any anomaly in Eq at the supposed phase 
transition point, nor would we expect this. |— ■ 

The most direct way of identifying any critical point {t/J)c is from poles of Dlog Fade approximantsll3. However 
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FIG. 4: The ground state energy for the square lattice (a) and for simple cubic lattice (b). The solid lines are different orders 
of integrated differential approximants to the dimer expansion series, while the points with errorbars are the results from Ising 
expansions. 
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FIG. 5: The staggered magnetizations for both localized spins (Mi) and conduction electrons (Ale) for the square lattice (a) 
and for sirrmle cubic lattice (b). Also shown for square lattice are the staggered moments m obtained from Monte Carlo 
calculationsQ. 
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FIG. 6: The inverse antiferromagnetic spin susceptibilities for both local (xi) and intinerant spins (xc) for the square lattice 
(a) and for simple cubic lattice (b). 



the series are irregular (Table I) and, perhaps not surprisingly, this yields no consistent results. However if we know, 
or assume, the value of the critical exponent 7 then biased estimates of the critical point can be obtained from direct 
Pade approximants to the series for x^/''', which should have a simple pole. Here we expect the transition for the square 
lattice to be of the same universality class as the d = 3 classical Heisenberg model with 7 ~ 1.4, and for the simple 
cubic lattice to be in the d = 4 universality class, with 7 = 1.0. In Table H we show estimates of x'^ {x = t/J) obtained 
in this way. As can be seen, a number of consistent estimates of the pole are obtained, particularly for the series for 
the conduction electron susceptibility. We might reasonably estimate for the square lattice [t/ J)c — 0.68 ± 0.02, and 
for the simple cubic lattice {t/J)c — 0.46 ± 0.01, where the error are subjective confidence limits. These correspond 
to {J/t)c = 1.48, 2.15 for the 2D and 3D case, values which are consistent with previous estimates. 

An alternative approach is to evaluate the staggered magnetization and susceptibility directly via integrated differ- 
ential approximants and to look at the behaviour as a function oft/ J. The staggered magnetizations shown in Figure 
H are obtained from an Ising expansions, starting from an antiferromagnetically ordered state. These are relatively 
constant for t/J > 1, but drop sharply to zero around t/J^ 0.6 — 0.7 in 2D and t/J^ 0.5 in 3D, confirming the 
existence of a transition to a magnetically disordered phase. The error limits are rather large, and it is not possible 
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FIG. 7: The triplet spin excitation spectrum (a) and quasiparticle excitat*^ spectrum (b) for the square lattice. The open 
points with errorbars are variational Monte Carlo results for a 6 x 6 latticeT 
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FIG. 8: The triplet spin gap As/ J and quasiparticle gap Aqp/J for the square lattice (a) and for simple cubic lattice (b) 
obtaii^d from different orders of integrated differential approximants. The points are the results of a Quantum Monte Carlo 
studyQ. 



to determine the transition point with high precision in this way. Our magnetization curves in 2D are very similar to 
the Quantum Monte Carlo resultsQ, although our conduction election magnetization is much smaller. We show the 
QMC results for comparison. In Figure |^ we show curves of the inverse susceptibilities, in the Kondo phase, obtained 
from dimer expansions. Again there is clear evidence for a transition, but it is difficult to locate precisely. Our best 
estimates from these figures would be t/J ~ 0.75 ± 0.10 in 2D and t/J ~ 0.50 ± 0.05 in 3D, less precise but quite 
consistent with the direct Fade approximant estimates. 

Next we consider the spin triplet and 1-hole (quasiparticle) excitations. Figure |^(a) shows the full spin triplet 
dispersion curve for t/ J = 0-25, 0.4, both in the Kondo phase, for 2D. For comparison we also show variational Monte 
Carlo results of Wang et alS3. The lowest spin excitation is at (tt, tt), and the spin gap clearly decreases as t increases. 

In Figure ^(b) we show the dispersion relation for the 1-hole excitation, again for t/J — 0.25, 0.4, for the 2D case. 
The minimum occurs at k = (0,0), with a maximum at (7r,7r). The overall shape is qualitatively similar to the 
mean-field results of Ref. We note the large error bars near (7r,7r) for t/J-. = 0.4. This may be a signature of 
degeneracy with a 2-particle continuum, as was observed in the Ji — J2 modeltj. 

Series have also been computed directly for the spin gap and quasiparticle gap, for both 2D and 3D cases. Analysis 
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is shown in Figure 8. The data clearly show the spin gap decreasing to zero at a critical point {t/J)c, the position 
being consistent with our estimates above. For the square lattice our results are in excellent agreement with the QMC 
results. For the 3D lattice we are unaware of any previous results for either energy gap. 



IV. SUMMARY AND DISCUSSION 



We have used linked cluster series methods in a comprehensive study of zero-temperature properties of the Kondo 
lattice model at half-filling, far the linear chain, square lattice, and simple cubic lattice. Our work significantly 
extends a previous series studyO for the ground state energy and susceptibility, and presents new series results for the 
magnetization, spin and quasiparticle dispersion relations, and energy gaps. Wherever possible we have compared 
our results with calculations by other methods and, in general, find excellent agreement. 

Our analysis supports the existence of a quantum critical point in the 2D and 3D cases, separating a Kondo spin- 
liquid phase (large J) from an ordered phase. Our estimates for the critical point are {t/J) = 0.68 ± 0.02, 0.46 ± 0.01 
in 2D and 3D respectively. These results are from direct Fade approximants to the series for x^^'^j ^^^^d hence are 
biassed by, but not particularly sensitive to, the choice of 7. The critical point estimate in 2D agrees very well with 
previous estimates, while in 3D we find a slightly higher value of (J/t)c, 2.17, compared with 1.833 from mean-field. 

We have not yet analysed the data for the frustrated close-packed triangular and face-centered cubic lattices, or for 
the body-centered cubic lattice, which we also have computed. Nor have we explored^te possible existence of bound 
states, or the region away from half-filling. Recent developments in series methodsEZrE3 make this possible, and we 
intend to pursue these directions, as well as others, in future work. 



Appendix 

The dimer series of ground state energy for ID KLM are: 

Eo/NJ = -3/4 - 2/3A2 - 14/45A'' + 0.7114685479^^^ + 0.19256921565A® 

-2.8528410569A^" + 2.2809484235A^2 + 12.882850521A^'^ - 30.446998097A^*^ 
-43.303667770A1® + 270.26259775A^*' + 0{\^'^) (9) 

The triplet spin excitation spectrum for ID KLM are 

8 A"^ 2356 A^ 

A^(fc) / J = 1 + 181.93078 A^ -f 2152.5066 A*^ - 27691.038 + 329428.5445 A^^ 

3 135 

/ 272 A^ 

-2.5917002 10^ A^'' - 2.75912 10^ A^^ H- 2.23129 10^ A^^ -t- ( 4 A^ + 305.675 A^ - 4076.18 A** 
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+55158.531287 A^" - 694363. A^^ + 6.4442 10^ A^'' + 2.20385 10^ A^^ - 3.67419 10^ A^^) cos(fc) 

(308 A'^ 
203.122 A^ + 3199.476135 A* - 50109.7 A^*' + 729779.9656 A^^ _ 8.92778 10*^ A^'' 
27 

-H5.42326 10^ A^^ + 1.72156 10^ A^*) cos(2 k) + ^71.0803 A^ - 1809.21 A^ + 35912.6434 
-640789. A^^ + 1.0026 10^ A^'^ - 1.22082 10^ A^^ + 4.1027 10^ A^^) cos(3 k) + (^521.4000489 A*^ 

-17575. Ai° + 416938.0779 A^^ - 8.34562 10^ A^"* + 1.40949 10^ A^^ - 1.77087 10^ A^^) cos(4 k) 

+ (4305.106 A^" - 178962. A^^ -I- 4.91952 10^ A^"* - 1.09405 10^ A^*^ + 2.00128 10^ A^^) cos(5 k) 
(38024.843 A^^ _ 1.87593 10^ A^'' + 5.85515 10^ A^^ - 1.43507 10^ A^^) cos(6 k) 

+ (351790.5722 A^"* - 2.00559 10^ A^^ 7.00211 10** A^*) cos(7 k) 

(3.36497 10*^ A^^ - 2.17486 10® A^*) cos(8 k) + 3.30089 10^ A^^ cos(9 k) (10) 
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TABLE I: Series coefficients for dimer expansions for the ground energy per site Eq/JN, the minimum triplet spin gap As/ J, 
and the antiferromagnetic spin susceptibihties for both local and conduction spins {xi and Xc)- Nonzero coefficients {t/J)" up 
to order n = 12 for square lattice and simple cubic lattice are listed. 



n 


Eo/JN 


As/ J 




XI 






square lattice 









-7.500000000x10"^ 


1.000000000 


5.000000000x10"^ 


5.000000000x10"^ 


2 


-1.333333333 


-1.333333333x10^ 


1.777777778 


6.518518519 


4 


-8.888888889x10"^ 


1.611851852x10^ 


4.632888889 


3.850916872x10^ 


6 


4.231487360 


-2.640795767x10^ 


6.236968731 


4.067531551x10^ 


8 


-1.519899530x10^ 


-1.072124632x10'^ 


4.340134129x10^ 


6.687413797x10^ 


10 


8.900723089x10"^ 


5.495699137x10'' 


1.004232541x10^ 


1.823289007x10^ 


12 


3.109971534x10^ 


-1.561259011x10* 


-6.384897515x10^ 


-4.654864554x10^ 






simple cubic lattice 









-7.500000000x10"^ 


1.000000000 


5.000000000x10"^ 


5.000000000x10"^ 


2 


-2.000000000 


-2.000000000x10^ 


2.666666667 


9.777777778 


4 


6.666666667x10"^ 


3.063111111x10^ 


1.291555556x10^ 


8.443753086x10^ 


6 


7.834807760 


1.060237771x10^ 


5.132587514x10^ 


2.399390035 xlO^ 


8 


-7.114791245x10^ 


-2.760658123x10^ 


3.510946982x10^ 


1.684588635x10^ 


10 


3.498368486x10^ 


3.860398011x10" 


1.354232378x10'^ 


1.217247537x10" 


12 


-1.603723348x10^ 


5.557190900x10* 


2.797156237x10^ 


-2.514875289x10* 
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TABLE II: Estimates of = {t/J)c from poles of [N,D] Pade approximants to the series for x, 
series for conduction electron, localized spin respectively. 



1/7 

c,i ■ 



The index c, I denotes the 



Square Lattice 


Approximant 




Approximant 




7 = 1.4 


[2,3]c 


0.45378 


[3,3]c 


0.46298 




[3,3]i 


0.46990 


[2,4]e 


0.46350 




Estimate xl 


= 0.46 ± 0.02, {t/J)c = 0.68 ± 0.02 
2 




2 



7 = 1.0 



[1,11c 

[2, 2], 
[l,3]c 
[2,3]c 
[3,3]c 
[2,4]. 
Estimate 



0.21 ±0.01 



0.20647 [0, 2]c 

0.21660 [2,2]i 

0.22094 [3, 2]c 

0.20601 [4, 21c 

0.20798 [3, 3]i 
0.20809 

, {t/J)c = 0.46 ±0.01 



0.20888 
0.21637 
0.19300 
0.19480 
0.20078 



